Analytic modelling of tidal effects in the relativistic inspiral of binary neutron stars 



O 

(N 

O 

Q 

(N 
(N 



o : 

i 

6fJ. 



> 
(N 

in 
O 

ON 

o 
o 



X 



Luca Baiotti, 1 ' 2 Thibault Damour, 3,4 Bruno Giacomazzo, 5 - 6 ' 7 Alessandro Nagar, 3 and Luciano Rezzolla 7 

1 Institute of Laser Engineering, Osaka University, Suita, Japan 
2 Yukawa Institute for Theoretical Physics, Kyoto University, Kyoto, Japan 
3 Institut des Hautes Etudes Scientifiques, Bures-sur-Yvette, France 
4 ICRANet, Pescara, Italy 
' 'Department of Astronomy, University of Maryland, College Park, Maryland USA 
b Gravitational Astrophysics Laboratory, NASA Goddard Space Flight Center, Greenbelt, Maryland USA 
1 Max-PIanck-Institut fur Gravitationsphysik, Albert-Einstein-Institut, Potsdam, Germany 
8 Department of Physics and Astronomy, Louisiana State University, Baton Rouge, Louisiana, USA 

To detect the gravitational-wave (GW) signal from binary neutron stars and extract information about the 
equation of state of matter at nuclear density, it is necessary to match the signal with a bank of accurate templates. 
We present the two longest (to date) general-relativistic simulations of equal-mass binary neutron stars with 
different compactnesses, C = 0.12 and C = 0.14, and compare them with a tidal extension of the effective-one- 
body (EOB) model. The typical numerical phasing errors over the ~ 22 GW cycles are A</> ~ ±0.24 rad. By 
calibrating only one parameter (representing a higher-order amplification of tidal effects), the EOB model can 
reproduce, within the numerical error, the two numerical waveforms essentially up to the merger. By contrast, 
the third post-Newtonian Taylor-T4 approximant with leading-order tidal corrections dephases with respect to 
the numerical waveforms by several radians. 
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Introduction. Inspiralling binary neutron stars (BNSs) are 
among the strongest sources of gravitational waves (GWs) and 
certain targets for the advanced and new-generation ground- 
based GW detectors LIGO/Virgo/GEO/ET H. These detec- 
tors will be sensitive to the inspiral GW signal up to GW fre- 
quencies of about 1 kHz, which are reached just before the 
merger. The late inspiral signal will be influenced by tidal 
interaction between the two neutron stars (NSs), which, in 
turn, encodes important information about the equation of 
state (EOS) of matter at nuclear densities. However, to re- 
liably extract such information, both a large sample of numer- 
ical simulations and an analytical model of inspiralling BNSs 
which is able to reproduce them accurately, are needed. In this 
Letter we report on significant progress on this problem by 
presenting the longest (to date) simulations of merging equal- 
mass BNSs and by showing how to use them to calibrate an 
effective-one-body (EOB) model of tidally interacting BNSs. 

Numerical simulations of merging BNSs in full general rel- 
ativity have reached a high-level of accuracy and have become 
more realistic (e.g. , with the inclusion of magnetic fields) 
only recently J3-Q11- The analytical description of tidally - 
interacting binary systems has been initiated very recently |a- 
0], with the following two major results. First, the dimension- 
less quantity kg (Love number) in the (gravito-electric) tidal 
polarizability parameter Gyn = 2kiR 2l+1 / (21 — 1)!! measur- 
ing the relativistic coupling (of multipolar order £) between a 
NS of radius R and the external gravitational field in which it 
is embedded has been found to be a strongly decreasing func- 
tion of the compactness parameter C = GMj (c 2 R) of the NS. 
Second, a comparison between the numerical computation of 
the bin ding energy of quasi-equilibrium circular sequences of 
BNSs and the EOB description of tidal effects H has 
suggested that higher-order post-Newtonian (PN) corrections 
to tidal effects effectively increase by a factor of order two 
the tidal polarizability of close NSs. The main aim of this 



paper is to extend the domain of applicability of the EOB 
method ifTHl . from the inspiralling binary black hole (BBH) 
case (for which it recently provided a very accurate analytic 
description lfl2[[l3ll ). to the yet unexplored case of inspiralling 
BNSs. To this aim we have performed two accurate and long- 
term BNS simulations covering approximately 20 — 22 GW 
cycles of late inspiral, and we will show that they can be re- 
produced accurately almost up to the merger by a new tidal 
extension of the EOB model. 

Tidal corrections in the EOB approach. We recall that 
the EOB formalism ifTTIl replaces the PN-expanded two-body 
dynamics by a resummed description with, in particular, a 



Hamiltonian of the form: Peob = Mc 2 y 1 + 2v(H e f( — 1), 
where M = Ma + Mb is the total mass and where 
is = Ma Mb / (Ma + Mb) 2 is the symmetric mass ratio. 
Here the "effective Hamiltonian" H c f[ is a simple function 
of the momenta and it incorporates the relativistic gravita- 
tional attraction mainly through the so-called "EOB radial 
potential" A(r). The structure of A(r) is remarkably sim- 
ple at 3PN: A 3PN (r) = 1 - 2u + 2vv? + a 4 vu A , where 
a 4 = 94/3 - (41/32)tt 2 , and u = GM/(c 2 r AB )- An excel- 
lent description of BBHs has been found to be given by lfl2tl 
A°(r) = Pg 1 [1 - 2u + 2vu 3 + a 4 vu 4 + a b uu b + a 6 isu 6 ] , 
where P T ™ denotes an (n, to) Pade approximant and where 
values of the coefficient a$ = —6.37, ag = +50 provide a 
very good agreement between EOB and numerical-relativity 
(NR) waveforms for BBHs ll 2ll (the results presented here are 
insensitive to this choice as long as and are chosen in a 
well defined range). Ref. H] suggested to include tidal effects 
as corrections both to the radial potential and to the waveform 
(and radiation reaction). The tidally corrected radial potential 
reads A(u) = A°(u) + A tidal (u), where 
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Here kJu 21+2 describes the leading-order (LO) tidal interac- 
tions. It is a function of the two masses, of the two compact- 
nesses Ca,b, and of the two (relativistic) Love numbers k] 
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The additional factor A^ ldal (w) in Eq. (Q]) represents the ef- 
fect of higher-order relativistic contributions to the tidal inter- 
actions: next-to-leading order (NLO), next-to-next-to-leading 
order (NNLO), etc. A number of different prescriptions can 
be considered for the correcting tidal factor ^L^^ 1 and these 
will be presented in a longer companion work lfl4ll . Here, we 
will limit ourselves to an ^-independent, "Taylor-expanded" 
expression ^ ldal (u) = 1 + ol\U + 6t2U 2 HI, where a n are 
pure numbers in the equal-mass case, but functions of Ma, 
Ca, and kf in the general case. The analytical value of the 
{£ = 2) 1PN coefficient ct\ has been reported in U (and re- 
cently confirmed in lfl5ll ). In the equal-mass case, it yields 
a,\ = 1.25. We will use this analytical value in the following 
and use our simulations to constrain the value of the 2PN co- 
efficient a.2- Similarly, one takes into account an £ = 2 tidal 
correction to the waveform and radiation reaction, as given 
at LO in Sec. V of H]. Additional coefficients parametrizing 
higher-order tidal relativistic contributions in the waveform 
and radiation reaction (such as 0i in Eq. (71) of JUl), were 
found to have a small effect UM and will be neglected here. 
In principle, tidal effects can also be accounted for via modi- 
fications of one of the non- re summed PN models, such as the 
Taylor- T4 one; see below for its comparison with the NR re- 
sults. 

In order to measure the influence of tidal effects, it is useful 
to consider the "phase acceleration" u> = du/dt = d 2 (j)/dt 2 , 
where <fi = (f>22 is the phase of either the curvature or of the 
metric I = in = 2 GWs. The function ui(u) is independent 
of the two "shift ambiguities" that affect the GW phase 4>{t), 
namely the shifts in time and phase, and thus a useful intrinsic 
measure of the quality of the waveform lfl7ll . However, here 
we use another diagnostic to measure the phase acceleration, 
namely the dimensionless function 



d(f) ujd4>/dt lo 2 
rflnw diu/dt lu 
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Numerical Simulations. They were performed with the 
Cactus-Carpet-Whisky IU8I1 codes and, in essence, we 
use the same gauges and numerical methods already discussed 
in JH. As initial data we use quasi-equilibrium irrotational bi- 
naries generated with the multi-domain spectral-method code 
LORENE, within a conformally-flat spacetime metric lfl9tl . 
The EOS of the initial data is the polytropic one p = K p T , 
where p, p, K = 123.6, and T = 2 are the pressure, rest- 
mass density, polytropic constant, and adiabatic index, re- 
spectively (in units where c = G = Mq = 1). The evo- 
lutions are instead performed with either a polytropic EOS 
or an "ideal-fluid" one, p = pe(T — 1), where e is the spe- 
cific internal energy; the differences in phasing introduced 
by the different EOSs are of order ±0.13 rad 11411 . Because 
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FIG. 1 . Comparison of the EOB curves for different choices of 
the effective tidal amplification factor j4^ ldal (it) = 1 + ai« + a.2U 2 ', 
with the corresponding NR ones (dashed lines with open circles) 
for the two binaries considered. The dotted line corresponds to the 
"tidal-free" ( or "point-mass") EOB, namely when ignoring tidal ef- 
fects. The figure also includes two Taylor-T4 models: tidal-free, and 
augmented by LO tidal effects. 



the stellar compactness represents the most important param- 
eter determining the size of tidal effects, we have considered 
two different (equal-mass) binaries having total Arnowitt- 
Deser-Misner (baryonic) mass of either 2.69 (2.89) Mq or 
3.00 (3.25) M©, thus with compactnesses C = Ca = Cb = 
0.12 or C = 0.14. The dominant (£ = 2) tidal parameters for 
the two compactnesses C — 0.12 (0.14) are found to be 10], 
respectively, k 2 = k£ = fcf = 0.00969 (0.07894), and there- 
fore K2 = 496.01 (183.81). Hereafter the two binaries will 
be referred to as M2 . 9C . 12 and M3 . 2C . 14, respectively. 
The number of refinement levels and their resolutions are the 
same as those in J2l, but the initial coordinate separation be- 
tween the stellar centers is 60 km, i.e. considerably larger than 
the one considered in 0]. This yields about 10 orbits be- 
fore merger, thus the two longest BNS waveforms produced 
to date. 

Discussion. We start our comparison between the NR re- 
sults and the EOB ones by showing in Fig.[T]the Q u diagnos- 
tics for various possible LO/NLO tidal models and for scaled 
GW frequencies Mlo < 0.06 [i.e. up to 3 (5) GW cycles be- 
fore merger for the M2 . 9C . 12 (M3 . 2C . 14) binary]. [We 
have estimated the NR Q w by fitting each NR phase evolu- 
tion on a (scaled) frequency interval / = [0.043, 0.057] with 
an analytical expression that reproduces at lower order the be- 
havior expected from the PN approximation, thereby filtering 
out the amplified numerical noise coming from the two time 
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FIG. 2. Comparison between NR and EOB phasing for the M2 . 9C . 12 (left panels) and M3 . 2C . 14 (right panels) binaries. The top panels 
show the real parts of the /122 waveforms, while the bottom panels show the corresponding phase differences. Note the excellent agreement 
almost up to the time of merger (vertical dashed and dot-dashed lines) and the very large errors when tidal effects are neglected (dotted line). 



derivatives in u = d 2 <fi/dt 2 (more details will be presented 
in ifTilo i. The first thing to note is that both the tidal-free 
EOB model (EOB TF , dotted line) and the EOB model includ- 
ing only LO tidal corrections (EOB LO , dot-dashed line) are 
clearly unable, both for the M2 . 9C . 12 (upper panel) and the 
M3 . 2C . 14 binaries (lower panel), to match the correspond- 
ing NR curves (dashed line with open circles). The dephas- 
ing accumulated over the frequency interval / A/(/) EOBNR = 
jj(<35 OB - <2^ R )dlnw, by the EOB LO model relative to 
the C = 0.12 (0.14) NR data is about 5.5 (2.0) rad. This 
is much larger than the NR phasing error, related to a finite- 
radius extraction and EOS dependence, estimated to be A<fi = 
±0.24 01. 

The inclusion of the NLO 1PN tidal effect (a x = 1.25 H) 
only slightly reduces these dephasings to about 4.9 (1.8) rad 
(EOB NLO curves in Fig. [TJ. This clearly indicates the need 
for NNLO (2PN and higher) tidal effects. We then found 
that choosing d'2 ~ 130, yields a good match between the 
Q u curves (solid line, EOB NNLO ) and the NR data (dashed 
line with open circles) for both binaries, with a corresponding 
phase difference A/0 EOBNR w 0.1 rad. The value a 2 « 130 
is probably only an effective description of higher-order rel- 
ativistic tidal effects. Moreover, the precise value of 0*2, or 
more generally of the amplification factor A^ ldal (u), is sen- 
sitive to the numerical truncation error. When considering 
resolution-extrapolated GWs 11411 . we found a smaller value 
of a.2, which is compatible with the estimate obtained using 
the binding energy of circular BNSs lioll . 

Figure Q] also reports the Q u diagnostics obtained when 
using two versions of the Taylor- T4 approximant: the tidal- 
free model (T4 TF , magenta, upper solid line), and the Taylor- 



T4 LO one (thick-dashed line). The latter model includes 
only the LO tidal effects fl, i.e. the LO tidal contribution 
a tida \x) oc k%x 5 to dxldt [where x = (Afw/2) 2 / 3 ; see H 
and Eqs. (86)-(88) of Note that the tidal-free Taylor-T4 
Qu: curve nearly coincides with the tidal-free EOB one, with a 
dephasing A/0 T4EOB = 0.013 rad. On the other hand, the I- 
integrated dephasings between the T4 LO description and the 
NR results are rather large, namely A/</> T4NR = 6.96 (2.53) 
rad for C = 0.12 (0.14). We have investigated whether a suit- 
able PN-amplification factor d tldal (x) = 1 + aj 4 x + aj 4 x 2 
of a tldal (x) might accurately reproduce the NR data, but 
we found that this was not possible with a single choice of 
d tldal (a;) for the two simulations Hill . The latter result sug- 
gests that the EOB modelling of tidal effects may be more 
robust than the corresponding Taylor- T4 one 

We next consider the comparison of the waveforms in the 
time domain and over the full inspiral up to the merger. This 
is shown in Fig. [2] whose left panels refer to the M2 . 9C . 12 
binary and the right ones to M3 . 2C . 14. The top parts com- 
pare the (real part) of the EOB NNLO (with a x = 1.25, a 2 = 
130) and NR metric /122 waveforms, while the bottom panels 
show the corresponding phase differences, A(/) EOBNR (i) = 
EOB (i) — NR (i) (suitably shifted in time and phase a 
la UH]). The two vertical lines indicate two possible markers 
of the "time of the merger"; more specifically, the dashed lines 
refer to the NR merger, defined as the time at which the mod- 
ulus of the metric NR waveform reaches its first maximum, 
while the vertical dash-dotted line represents the EOB esti- 
mate of the "formal" contact U. Figure [2] clearly shows that 
the agreement in the time domain between the analytic EOB 
description and the numerical one is extremely good essen- 
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tially up to the merger, with a phase error which is well within 
the estimated error bar. More precisely: (i) in the M3 . 2C . 1 4 
case, the phase difference At/> EOBNR (i) varies between —0.1 
rad in the early inspiral and +0.04 rad at the NR merger; (ii) in 
the M3 . 2C . 12 case, A(/) EOBNR (i) varies between -0.15 rad 
in the early inspiral and —0.15 rad as late as 130M (i.e. ap- 
proximately 1.5 GW cycles) before the NR merger. For the 
latter binary, it is only just before the NR merger that the de- 
phasing becomes of order 1 rad. 

Conclusions. We have presented the first NR-EOB com- 
parison of the GWs emitted during the inspiral of relativis- 
tic BNSs. In particular, we have analyzed the longest to date 
numerical simulations of equal-mass, irrotational BNSs with 
two different compactnesses. We found that tidal effects are 
significantly amplified by higher-order relativistic corrections 
(2PN and higher) and that the inclusion of such corrections is 
necessary for an accurate phasing of the GW signal. Such an 
amplification is equivalent to a (distance-dependent) effective 
increase of the Love numbers. When a single choice for the 
unique free parameter in the NNLO term is made, the EOB 
model is able to reproduce the two NR phase evolutions well 
within the estimated NR error and essentially up to merger 



(and in fact up to merger in the C = 0.14 case). By contrast, 
we have shown that the use of the Taylor- T4 PN approximant 
considered in ^ leads to phase disagreements (over the fre- 
quency interval / = [0.043, 0.057]) A 7 </> T4NR = 6.96 (2.53) 
rad for C = 0.12 (0.14). 

The work reported here provides the first evidence that an 
accurate analytic modelling of the late inspiral of tidally in- 
teracting BNSs is possible, thereby opening the possibility to 
extract reliable information on the EOS of matter at nuclear 
densities from the data of the forthcoming advanced GW de- 
tectors. These encouraging results, however, also call for a 
continued synergy between more accurate numerical simu- 
lations (notably exploring different mass ratios) and higher- 
order analytic results. 
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